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摘 要 针对 水 平 管道 中 油水 分 散 流 流速 测量 问题 , 提出 一 种 基于 超声 互相 关 的 测量 方法 . 通过 在 水 平 管道 上 间隔 一 定 
距离 的 上 下 游 安装 两 对 超声 探头 获取 油水 两 相 流 流动 过 程 中 的 随机 波动 信号 , 针对 离散 相 和 连续 相 之 间 的 相对 滑动 , 将 互 
相关 方法 与 漂移 模型 和 运动 波 理 论 相 结 合 , 并 考虑 流速 分 布 和 含 率 分 布 对 流速 测量 带 来 的 影响 , 在 连续 相 不 同 的 情况 下 分 
别 建立 两 相 流 总 表 观 流速 的 测量 模型 , 并 采用 油水 两 相 流 实际 流动 状态 的 实验 数据 计算 模型 参数 . 结果 表明 , 测量 模型 对 
流 型 具有 依赖 性 , 漂移 速度 是 离散 相 含 率 的 函数 . 通过 实验 数据 对 模型 进行 验证 , 针对 不 同 流 型 选取 不 同 参数 计算 得 到 的 
总 表 观 流速 与 实际 流速 能 够 较 好 吻合 , 相对 误差 在 6.12% 以 内 . 
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Abstract For the overall superficial velocity measurement of dispersed oil-water flow in horizontal 
pipeline, a measuring method based on ultrasonic cross-correlation technique is introduced. Two 
pairs of ultrasonic probes are installed in the upstream and downstream within a certain distance to 
obtain the flow signal of two-phase flow. Considering the slip between the dispersed phase and the 
continuous phase and the influence of velocity profile and phase distribution on the measurement, 
the overall superficial velocity measurement model is established by combining the cross-correlation 
algorithm, drift-Hux model and kinematic wave theory. The parameters of this model were calculated 
by the experimental data under different flow states. The results show that the measurement model 
= is dependent on the flow regime, and the drift velocity is a function of the holdup. The model is 
validated by experimental data, and it shows that the overall superficial low velocity calculated by 
this model is in good agreement with the actual How velocity. The relative error is less than 6.12%. 
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0 引言 


油水 两 相 流 广 泛 存在 于 能 源 、 化 工 等 工业 生产 目前 对 于 油水 两 相 流 流速 测量 的 研究 已 有 很 多 . 
以 及 石油 的 开采 、 运 输 、 存 储 过 程 中 ,根据 连续 相 。 互相 关 方 法 是 一 种 常用 的 流速 测量 方法 ， 它 通过 测 
的 不 同 ， 可 将 其 分 为 油 连 续 和 水 连续 两 种 形式 趾 。 ” 量 流 动 噪 声 在 从 一 个 传感器 传 到 另 一 个 传感器 时 所 
油水 两 相 流 流速 的 准确 测量 对 预 佑 产能 、 保 障 生产 。 用 的 渡 越 时 间 ， 在 传感器 间距 已 知 的 情况 下 进而 得 
安全 以 及 研究 管道 中 流体 的 流动 特性 具有 重要 的 意 。 到 流体 流速 ， 可 以 使 用 光学 、 电 学 、 声 学 等 传 感 
义 。 但 是 由 于 油水 两 相 不 同 的 介质 属性 ， 两 相 流 的 咒 同 。Lucas 等 是 利用 电阻 互相 关 流 量 计 研究 了 
流动 过 程 呈现 了 复杂 的 、 高 随机 性 的 流动 状态 , 流 ， 倾斜 管道 下 油水 两 相 流 互 相关 流速 与 平均 流速 、 倾 
型 变化 尤为 复杂 ， 且 油 水 相间 存在 着 相对 滑动 中 ， ” 角 及 含 率 的 关系 。Tan 等 『- 引 分 别 利 用 单 层 、 双 层 
这 些 都 为 实际 测量 带 来 了 困难 和 挑战 B- 委 。 电阻 层 析 成 像 系 统 和 互相 关 技 术 对 两 相 流 流 速 进 行 
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3 期 史 雪 芽 等 : 油水 分 散 流 超声 互相 关 流 速 测量 方法 


了 测量 . Zhai 等 四 利用 双 电 容 探 针 测量 了 水 平 管 油 
水 两 相 流 的 互相 关 流 速 . Zhou 等 9 将 电容 耦合 式 
非 接触 电导 测 基 结 合 互 相关 方法 用 于 毫米 级 管 径 气 
液 两 相 流 的 气泡 速度 测量 。 相对 于 电学 传感器 而 言 ， 
由 于 超声 波 在 传播 过 程 中 对 流 场 不 产生 扰动 、 无 压 
力 损 失 ， 可 实现 非 侵入 的 参数 检测 ， 并 且 其 安全 性 
更 有 保障 , 是 流动 参数 检测 的 一 种 重要 手段 9， 

两 相 流 流速 测 其 多 采用 均 相模 型 建 模 , 但 在 实 
际 流动 过 程 中 ,两 相 流 相间 存在 着 随机 可 变 的 相 界 
面 和 相间 滑动 ， 均 相模 型 并 不 能 满足 测量 精度 要 
求 。Zuperba 在 1965 年 提出 漂移 模型 的 概念 ， 阅 
明了 相间 相对 速度 和 总 表 观 流速 、 相 含 率 之 间 联 系 ， 
模型 综合 考虑 了 流体 的 不 均匀 分 布 和 相间 滑动 对 流 
速 测量 带 来 的 影响 。 模型 中 有 两 个 重要 参数 : 分布 
参数 和 漂移 速度 ， 前 人 对 此 做 了 很 多 研究 。Zuber 
等 9 得 出 相 含 率 大 小 会 影响 分 布 参 数 的 取 值 ， 漂 
移 速度 受 流 型 影响 。Ishii 等 [04 认为 分 布 参数 与 密 
度 比 和 雷诺 数 有 关 ， 并 利用 流体 的 运动 方程 和 连续 
性 方程 , 考虑 几何 体 、 前 应力、 界面 的 动量 传递 等 因 
素 的 影响 ， 在 弹 状 流 、 搅 拌 流 和 环 状 流下 的 漂移 速 
度 进行 了 计算 . Hibiki 等 3 考虑 气泡 尺寸 的 影响 ， 
研究 了 泡 状 流下 的 漂移 模型 . Lucas 等 9] 考虑 管道 
倾斜 角度 的 影响 ， 对 模型 参数 进行 了 讨论 . 可 以 看 
出 ,漂移 模型 的 参数 受 表面 张力 、 雷 诺 数 、 剪 应力 、 
倾角 、 含 率 分 布 、 流 型 等 因素 影响 , 模型 参数 的 确定 
对 流速 的 准确 计算 具有 重要 意义 。 

因此 , 针对 水 平 管道 油水 分 散 流 流速 测量 问题 ， 
提出 一 种 基于 超声 互相 关 的 流速 测量 方法 。 考虑 油 
水 相间 的 滑动 现象 ,在 漂移 模型 的 基础 上 结合 运动 
波 理论 ,综合 考虑 流 型 和 含 率 分 布 对 流速 测量 带 来 
的 影响 ， 在 连续 相 不 同 的 情况 下 分 别 建立 两 相 流 总 
表 观 流速 的 测 基 模型 。 其 中 ， 基 于 油水 两 相 流 的 实 
际 流 动 状态 对 分 布 参数 进行 了 理论 推导 ， 漂 移 速度 
根据 实验 数据 获得 . 并 利用 实验 对 模型 进行 了 验证 ， 
提高 了 总 表 观 流速 的 测量 精度 ， 


1 互相 关 测 量 原理 


超声 波 互 相关 原理 如 图 1 所 示 。 两 对 超声 探头 
间隔 50 mm 分 别 安装 在 管道 的 上 下 游 , 其 中 51, 52 
为 超声 波 发 射 探 尖 ，R1，R2 为 超声 波 接收 探头 . 由 
于 油水 两 相 流 的 流动 是 一 个 随机 过 程 ， 超 声波 流 过 
管道 时 ， 声 波 信 号 会 受到 流体 速度 场 的 调制 。 根据 
“凝固 ”流动 图 形 假设 ,， 当 流体 中 存在 某 一 扰动 的 情 
况 下 , 经 过 调制 的 流动 信号 将 在 上 下 游 重复 出 现 , 并 
引起 相似 的 作用 , 只 是 时 间 上 滞后 . 因此 对 上 下 游 信 


号 进行 互相 关 则 可 以 计算 出 流速 。 上 下 游 信号 的 互 
相关 为 : 


T 
Ry (7) = 去 Z (ty (十 7T)d7 (1) 


式 中 , z(t), y(t) 为 上 下 游 的 随机 波动 信号 , 了 为 测 
量 时 间 。 则 随机 噪声 信和 号 传播 速度 可 表示 为 : 

Wec = L/To (2) 
式 中 , 工 为 两 对 传 感 需 之 间 的 距离 ，m 为 互相 关 函 
数 的 峰值 点 所 对 应 的 时 间 ， 即 为 流体 从 上 游 传感器 
到 下 游 传 感 器 所 用 的 湾 越 时 间 。vwec 为 互相 关 速度 ， 
通常 被 认为 是 混合 流体 的 结构 传播 速度 。 


发 射 传感器 


流动 方向 
一 


接收 传感器 


RI 


Z/m 


图 1 互相 关 原 理 


Fig. 1 _ Principle of cross-correlation method 


2 两 相 流 模型 
2.1 漂移 模型 和 运动 波 理论 

油水 两 相 流 中 不 同 的 相 之 间 存 在 着 相对 滑动 ， 
因此 采用 Zuber 和 Findlay 在 1965 年 提出 的 漂移 
模型 H3. 模型 形式 如 下 : 

Va = Co (7) + Vay (3) 

式 中 , Ta 为 离散 相 加 权 平 均 流 速 , (7) 截面 平均 的 混 
合流 体 的 表 观 流速 . Co 为 分 布 参数 , 5a; 称 为 加 权 平 
均 漂 移 速 度 . 分 布 参数 可 利用 下 式 进行 计算 : 


1 
A 


= (4) 
四 人 Th 
fl fm 
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式 中 , 4 为 管道 截面 面积 , a 为 离散 相 的 相 含 率 ，vd 
为 离散 相 漂 移 速 度 ， 在 一 般 的 工程 实际 条 件 下 ， 可 
认为 zu = va 站。 加 权 平 均 漂移 速度 ru 受 多 种 因素 
影响 , Zuber 指出 va; 满足 : 


Da; = Vo (1 — oo (5) 


式 中 , vo 为 单个 液 滴 的 终端 流速 ,k 一 般 取 2. 若 不 
考虑 各 相 的 变化 时 ,Zuber 和 Findlay 根据 流体 的 连 
续 性 方程 和 运动 波 理论 ， 指 出 运动 波 的 传播 速度 和 
漂移 模型 具有 一 样 的 函数 形式 ， 并 指出 由 互相 关 法 
获得 的 互相 关 速 度 wcc 就 是 运动 波 的 传播 速度 ux， 
得 到 : 

d (Co (7) + va;) 
Uec = Co 7 太 十 al) (7) 


式 中 , f(a) 称 为 运动 波 漂移 速度 . 因此 在 获得 Co 和 
f(a) 的 情况 下 ,可 通过 互相 关 流 速 得 两 相 流 体 的 总 
表 观 流速 。 
2.2 分 布 参数 的 确定 

水 平 管 油水 两 相 流 流 型 多 变 ， 并 且 在 含水 率 不 
同 的 情况 下 连续 相 也 不 同 。 对 于 水 连续 的 油水 分 散 
流 , 其 主要 流 型 包括 水 包 油 (o/w)、 油 包 水 和 水 包 油 
(Dw/o&Do/w)、 水 包 油 和 水 (Do/w&ew)。 而 对 于 油 
连续 的 油水 分 散 流 ,实验 中 能 观察 到 的 流 型 有 油 包 
水 (w/o)、 油 包 水 和 水 包 油 (Dw/o&Do/w)。 在 连续 
相 不 同 的 情况 下 , 流速 分 布 满足 不 同 的 关系 式 : 


Wec = Upw = Co (7) Va 二 a 


于 
H nn 

2aoR 十 Q2 (2 ) 
C0 = 


由 式 (10)、 式 (11) 可 得 , 当 五 =R 时 , 00 = 1; 


1 


水 连续 : 二 一 (1 一下)” (8) 
滑 连 续 , 二 -一 1 一 ( 主 ) (9) 


式 中 , j 为 总 表 观 流速 , jwax 管道 中 心 最 大 速度 ,7 为 
管道 任意 一 点 到 管道 中 心 的 距离 ，R 为 管道 半径 ,mw 
一 般 取 值 为 6~10， 当 n=7 为 符合 大 多 数 清流 状态 
下 的 经 验 取 值 。 

由 于 分 布 参数 受 流速 分 布 和 含 率 分 布 的 影响 ， 
针对 油水 分 散 流 ,根据 式 (4) 分 别 计算 在 不 同 的 流 
型 下 的 分 布 参数 值 . 

在 均匀 的 o/w 和 均匀 的 w/o 情况 下 ,由 于 离散 
相 的 分 散 程度 较 高 ， 可 以 忽略 离散 相 含 率 分 布 对 计 
算 带 来 的 影响 , 认为 Co=1。 

在 Do/w&w 流 型 下 , 油 相 和 水 相 之 间 存 在 着 相 
分 界面 . 设 油 水 两 相 流 分 层 的 层 高 为 五 (下 文 将 分 为 
五 > R 和 卫 < RR 两 种 情况 进行 讨论 )， 上 层 中 离散 
相 的 含 率 为 oil, 下 层 离散 相 的 含 率 为 a2。 根据 流速 
分 布 和 式 (4), 并 将 面积 分 简化 为 线 积分 , 在 五 > RR 
与 了 < RR 两 种 情况 下 分 别 得 到 Do/w&w 流 型 时 分 
布 参数 的 表达 式 : 

当 万 < 时 : 


本 
HN\r7 
(oa2 Pe ai)H (#) 十 2Rai 


Co = 10 
9 a2H + 2Rat — Has Wo 
当 万 > 忆 时 : 
(H -2R) + oR (2R— H)!t* 
11 
Q2H+2Ra 一 万 ai ( ) 
EE 
、 、 2R\™ 
六 连续 + ov = ( 到 (4) 


设 al = 1，a2 = 0, 分 别 将 两 式 化 简 ， 当 H/R 取 
0.1~2 时 分 别 计算 分 布 参数 的 值 ， 并 由 几何 关系 确 
定 有 /RR 和 含油 率 的 关系 , 可 以 得 到 分 布 参数 与 含油 
率 的 关系 Co = f(aQ) 为 : 


Co = —0.5120? + 0.758a + 0.757 (12) 


Dw/o&Do/w 流 型 在 水 连续 和 油 连续 时 均 会 出 
。 根据 Elseth 的 研究 071， 在 连续 相 不 同 的 情况 
， 离 散 相 径 向 含 率 分 布 分 别 满足 以 下 关系 式 : 


水 连续 ,oo = ( 汉 )” (13) 


才 党 


式 中 ,ao 为 含油 率 ，auw 为 含水 率 , y 为 流 道 任意 一 
点 到 管 底 部 之 间 的 距离 ,m 一 般 取 值 1~10. 

根据 经 验 取 值 ，m 取 3~8 时 根据 式 (4) 分 别 计 
算 连 续 相 不 同时 分 布 参 数值 , 结果 见 表 1, 表 2. 


表 1 水 连续 下 (”= 7) 分 布 参 数 取 值 表 
Table 1 The value of Distribution parameter 
(water continuous and n= 7) 

m 3 4 5 6 7 8 
Co 1.008 1.008 1.007 1.006 1.005 1.005 


PE 


a 


3 期 


表 2 油 连 续 下 分 布 参数 取 值 表 


Table 2 The value of Distribution parameter (oil 


continuous) 
mm 3 4 5 6 7 8 
Co 0.900 0.935 0.952 0.963 0.969 0.974 


2.3 运动 波 漂 移 速 度 的 确定 

运动 波 漂移 速度 f(a) 受 多 种 因素 影响 ， 确 定 
较为 复杂 。 根据 式 (5) 和 对 分 布 参数 的 理论 推导 可 
知 ,f (a) 是 离散 相 含 率 a 的 二 次 函数 ,可 表示 为 : 


f(a)=aao*+batc (15) 


式 中 , a, 0, c 为 待定 参数 ， 

因此 本 文 利用 实验 数据 采用 拟 合 的 方法 确定 运 
动 波 漂移 速度 f(a). 根据 式 (7) 可 知 , 具体 思路 是 : 
在 连续 相 不 同 的 情况 下 ,计算 不 同 流 型 时 的 分 布 参 
数 , 然后 在 超声 实验 数据 的 基础 上 , 计算 出 互相 关 流 
速 ucce, 并 认为 (7) 为 实验 中 两 相 流 总 表 观 流速 ， 则 
可 根据 下 式 来 计算 不 同 流 型 下 f (a) 的 取 值 : 


Ja) = vec 一 Co (7) (16) 


根据 式 (15), 将 f(a) 与 离散 相仿 率 a 进行 二 
次 拟 合 ,确定 待定 参数 a, 5, c 的 取 值 ， 则 可 以 得 出 
不 同 流 型 下 f (a) 的 表达 式 。 


3 实验 研究 

油水 分 散 流 的 流速 测量 实验 于 天 津 大 学 多 相 流 
实验 室 完成 , 实验 装置 结构 如 图 2 所 示 . 实验 用 水 为 
自来水 (密度 998 kg/m3, 动力 严 度 0.001 Pa.s), 油 
为 工业 白 油 (密度 800 kg/m3, 动力 黏度 0.017 Pa-s)。 
水 平 管 道 由 不 锈 钢管 组 成 ， 内 径 为 50 mm。 入口 与 
出 口 之 间 的 距离 为 16.6 m。 在 入口 处 设置 单 相 流量 
计 ， 其 测量 精度 为 土 0.5%。 组 合 传 感 器 安装 于 入 口 
处 下 游 约 12 m 处 , 便于 使 流 型 充分 发 展 . 实验 中 固 
定 水 相 流量 ， 调 节 油 相 流量 ， 通 过 不 同 的 流量 配 比 
产生 不 同 流 型 , 其 中 油水 分 散 流 为 本 文 的 研究 对 象 . 
实验 中 两 相 流 总 流量 在 5.156 m3/h 到 17.960 m3/h 
内 变化 . 对 于 每 种 流动 状态 , 测量 过 程 持续 20 s。 

由 于 油水 两 相 流 的 流动 状态 较为 平稳 ， 因 此 可 
通过 快 关 阔 方法 获取 管道 内 流体 的 真实 相仿 率 。 快 
关 立 结构 如 图 3 所 示 。 在 每 一 流动 状态 下 ,同时 关 
闭 阀 门 1 和 阀门 2 对 流体 进行 截取 ,被 截取 的 流体 
经 过 静 置 后 ， 测 量 两 相 流 的 相 含 率 . 为 确保 测量 精 
度 和 排除 随机 误差 ， 实 验 中 每 一 流动 状态 下 对 流体 
截取 三 次 ， 取 三 次 测量 结果 的 平均 值 作为 该 流动 状 
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互相 关 流 速 测试 系统 结构 图 如 图 4 所 示 ， 两 对 
超声 探头 分 别 安装 于 管道 的 上 下 游 ， 上 下 游 传感器 
间隔 距离 为 50 mm, 信号 发 生 器 产生 1 MHz 连续 的 
激励 信号 , 经 过 功率 放大 后 送 入 超声 探头 的 发 送 端 ， 
经 过 流 场 的 调制 后 由 接收 端 进 行 接收 。 实验 中 采用 
PXI 信号 采集 板 卡 ， 结 合 LabVIEW 程序 设计 一 套 
数据 采集 上 位 机 软件 , 进行 数据 采集 , 采样 频率 为 5 
MHz, 采样 时 间 为 20 s. 对 获取 的 上 下 游 信号 进行 互 
相关 处 理 , 可 得 互相 关 流 速 。 
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图 2 实验 装置 结构 图 
Fig. 2 Flow facility 


图 3 快 关 疼 
Fig. 3 Quick closing values 


图 4 测试 系统 结构 图 
Fig. 4 System structure 


4 实验 结果 和 分 析 

将 两 相 流 总 表 观 流速 、 所 得 的 互相 关 流 速 以 及 
不 同 流 型 下 的 分 布 参数 值 代入 式 (16) 中 , 得 到 不 同 
流 型 下 的 f(a) 的 值 . 利用 式 (15) 在 不 同 流 型 下 将 
f(a) 与 离散 相 含 率 a 进行 二 次 拟 合 , 得 到 f(a) 与 
离散 相 含 率 a 的 关系 如 图 5、 图 6 所 示 。 需 指出 在 
Dw/o&Do/w 流 型 下 , 对 m 分 别 取 3~8 时 分 别 将 运 
动 波 漂 移 速度 f(a) 与 离散 相 含 率 a 进行 拟 合 ， 对 
比 发 现 m = 3 时 拟 合 程度 最 好 ,， 故 采用 此 结果 . 
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图 5 水 连续 下 f (a) 与 离散 相 含 率 a 的 关系 
Fig. 5 Relationship between thef (a) and the oil fraction 
Qwhen water is continuous 
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利用 实验 数据 ， 对 模型 进行 验证 ,将 不 同 流 型 
下 分 布 参数 Co、 运 动 波 漂移 速度 f(a) 以 及 互相 关 
流速 ucc 代入 式 (7), 计算 流体 的 总 表 观 流速 , 并 将 
模型 计算 出 的 总 表 管 流速 与 实验 测 基 的 流速 进行 对 
比 , 不 同 流 型 下 的 流速 测试 结果 如 图 7 所 示 ， 相 对 
误差 如 图 8 所 示 。 发 现 模 型 计算 结果 和 实验 测量 结 
果 吻 合 较 好 , 相对 误差 在 6.12% 以 内 , 平均 相对 误差 
为 1.15%, 证 明了 理论 模型 的 适用 性 和 准确 性 . 
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图 6 油 连 续 下 f (a) 与 离散 相 含 率 a 的 关系 
Fig. 6 Relationship between thef (a) and the oil fraction 
Qawhen oil is continuous 
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图 7 不 同 流 型 下 的 流速 测试 结果 


Fig. 7 Flow velocity test results under different flow regimes 


但 是 ,结果 中 仍然 存在 一 定 误 差 ， 来 源 主要 包 
括 以 下 几 个 方面 : 一 方面 ， 模 型 中 分 布 参数 的 计算 
依赖 于 对 流体 流速 分 布 和 含 率 分 布 的 假设 ， 但 是 由 
于 实际 流体 流动 复杂 ， 其 流速 分 布 和 含 率 分 布 并 不 
一 定 完全 符合 假设 , 并 且 n=7 只 是 能 符合 大 多 数 涪 
流 状 态 下 的 经 验 取 值 ， 男 一 方面 ， 互 相关 速度 的 实 
际 意义 并 不 明确 , 在 计算 时 的 计算 依赖 于 相关 时 间 ， 


3 期 


相关 次 数 的 选取 , 在 实际 计算 时 会 产生 一 定 误差 ; 此 
外 , 计算 中 的 一 些 简 化 运算 , 也 会 带 来 计算 方面 的 误 
差 . 因此 , 对 于 互相 关 速 度 的 物理 意义 、 流速 和 含 率 
的 分 布 仍 是 多 相 流 测量 中 所 要 解决 的 问题 。 
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Fig. 8 Relative error distribution under different flow regimes 


5 结 论 


针对 水 平 管道 中 油水 分 散 流 流速 测量 问题 ， 提 
出 一 种 基于 超声 互相 关 的 测 基 方法 . 通过 在 水 平 管 
道上 间隔 一 定 距 离 的 上 下 游 安装 两 对 超声 探头 来 获 
取 油 水 两 相 流 流动 过 程 中 的 随机 波动 信号 ,， 采用 互 
相关 方法 计算 两 相 流 互相 关 速 度 。 考虑 离散 相 和 连 
续 相 之 间 存 在 的 相对 滑动 现象 ， 将 互相 关 方 法 与 漂 
移 模型 、 运 动 波 理论 结合 ， 在 连续 相 不 同 的 情况 下 
分 别 建立 两 相 流 平均 混合 表 观 流速 的 测 基 模型 ， 在 
建 模 的 过 程 中 考虑 流 型 和 含 率 分 布 对 流速 测 基 带 来 
9 影响， 对 漂移 模型 中 的 分 布 参数 进行 了 理论 推导 ， 
并 通过 实验 数据 计算 了 漂移 速度 . 实验 表明 ,测量 
模型 对 流 型 具有 依赖 性 ， 且 漂移 速度 与 离散 相 含 率 
有 关 。 通 过 实验 数据 对 模型 进行 验证 ， 最 终 相 对 误 
差 在 6.12% 以 内 。 
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